The Yrast Spectra of Weakly Interacting Bose-Einstein Condensates 
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The low energy quantal spectrum is considered as a function of the total angular momentum for a 
system of weakly interacting bosonic atoms held together by an external isotropic harmonic potential. 
It is found that besides the usual condensation into the lowest state of the oscillator, the system 
exhibits two additional kinds of condensate and associated thermodynamic phase transitions. These 
new phenomena are derived from the degrees of freedom of "partition space" which describes the 
multitude of different ways in which the angular momentum can be distributed among the atoms 
while remaining all the time in the lowest state of the oscillator. 
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The properties of the non-zero angular momentum 
states of Bose-condensates of atoms in external traps 
have been addressed in several recent publications JlJ. 
Special interest attaches to the lowest energy quantum 
states with a given angular momentum. Borrowing from 
nuclear physics terminology we shall refer to this part 
of the quantum spectrum as the yrast region [2]. The- 
oretical work on the yrast region for Bose condensates 
has, with the exception of ||, been based on the Gross- 
Pitaevski mean field approximation [Q. In the present 
Letter the yrast spectra are described in terms of the el- 
ementary modes that carry the angular momentum. It is 
found that this approach reveals structures in the yrast 
spectra that were not recognized in previous work. The 
picture of the yrast region described in this note draws 
heavily on ideas developed in collaboration with A. Bohr 
H in connection with the analysis of the stability of per- 
sistent currents in superconductors. 

The model considered, which has previously been dis- 
cussed by Wilkin et al. H, involves N spin zero bosons 
moving in an isotropic harmonic confining potential and 
a contact interaction with a strength that is assumed to 
be weak, in the sense: 



vN < hu> , 



(1) 



where v is the expectation value of the interaction be- 
tween two particles both in the ground state of the oscil- 
lator, and huj is the quantum energy of the confining po- 
tential. Although the so-far published experiments have 
mainly studied systems that do not satisfy the condition 
(1), systems that do satisfy this condition are accessible 
with current experimental techniques and would be es- 
pecially interesting since this implies that the coherence 
length in the atomic cloud becomes larger than the size of 
the system, bringing into sharpest focus the mesoscopic 
nature of these systems. The eigenstates of a single par- 
ticle moving in an harmonic oscillator can be labeled by 
the quantum numbers (nlm) where n is the radial quan- 
tum number, I is the angular momentum, and m is the 



component of I on the axis of quantization. The energy 
of such a state is e = (2n + I + 3/2)Tlujq, and for the 
yrast states we need only consider occupancy of states 
with n — m = I > 0. (Without loss of generality we 
assume that the total angular momentum is in a state of 
maximum alignment L = M). Thus the relevant single 
particle states involve only a single non-negative quan- 
tum number, m. 

In the yrast state with total angular momentum TiL, 
the motion in the oscillator will contribute an excitation 
energy LKloq with respect to the ground state. How- 
ever, this configuration is degenerate in the absence of 
the interactions, since the oscillator energy is the same 
whether the L quanta are distributed among the parti- 
cles with one quantum to each of L different particles, or 
with 2 quanta on one particle and one quantum on each 
of L — 2 other particles etc. Since the multiplicity p(L) 
of the partitions of the integer L into positive integers 
increases asymptotically as M 



P(L) 



1 



L 



exp 




(2) 



the magnitude of the degeneracies becomes quite large, 
even for moderate values of L. This set of states will be 
called the partition space of L, and the problem of char- 
acterizing the yrast spectra thus amounts to finding the 
lowest quantum states selected by the contact interaction 
acting within the partition space of L. 

Partition space 

Since the construction of appropriate basis states in par- 
titions space and their matrix elements play the crucial 
role in the subsequent discussion, it is useful to begin 
with a brief discussion of some of the most elementary 
features of this interesting space. The standard notation 
for a partition of the integer L is 



(l ni 2 n2 ...L nL ), 



(3) 
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where the n r are restricted by the condition 

L 

J2 rn r=L. (4) 

r=l 

The set of integers {n r } satisfying (4) can thus be 
thought of as a set of quantum numbers characterizing a 
basis state in the partition space of L. 
In the present discussion we shall exploit two different 
ways of associating a basis state with a given partition 
{n r }. The two realizations differ in their identification of 
the object that carries the r units of angular momentum: 

(i) r = m, the address of a single particle state 

In this realization the n m count the number of different 
particles, each of which is carrying m units of angular 
momentum. 

*{« m } = \{m = l) ni (m = 2)" 2 .. .(m = L) nL ). (5) 

(ii) r = A, the multipole order of a (collective) normal 
mode of excitation 

Assuming that the normal mode is excited by acting on 
the ground state with the collective operator Q\ , we have 

* W} ~n£ =1 (Q A r lo>. (6) 

Although these two realizations each provide a com- 
plete set of states for describing partition space, these 
states can be very different when the numbers n r of ex- 
citations is large; thus for a single mode we have the 
overlap 




where we have employed the expression (8) for Q\. Thus 
for states involving condensates in which n\ <~ O(L) 
where L ^> \/]V, the two different states derived from 
the same partition become essentially orthogonal. 

In thinking about the structure of the operators Q\ it 
is fruitful to take inspiration from Feynman's arguments 
[6] that identify the long wavelength, low energy excita- 
tions of superfluid helium with modes that are excited by 
acting on the ground state with a sum of single particle 
operators that impart the conserved quantum numbers 
(momentum in the case of liquid helium, angular mo- 
mentum in the present problem) to a single particle of 
the condensate. Thus we shall assume 

1 N 

Q\ = ^=Y*t (8) 
V2AA!^ P 

Zp Xp -\- zypj (9) 



where the spatial coordinates are measured in units of 
the oscillator length. The operator (8) acting on the 
L = many body ground state of the oscillator creates a 
symmetrized state with angular momentum L = A; the 
operator (8) raised to the n\ power creates the many 
body excitations that appear in (6). 

The excitation energies of the normal modes (8) can 
be evaluated by taking expectation values of the contact 
interaction in the appropriate harmonic oscillator states 

e A = (OIO^FQaIO) - (0|V|0) (10) 

where the first term in the parenthesis represents the loss 
of interaction energy due to the removal of a single par- 
ticle from the condensate in the ground state (m = 0) of 
the oscillator, and the second term describes the remain- 
ing interaction between the condensate particles and the 
excited single particle moving in an orbit with angular 
momentum A. The short range of the interaction is re- 
sponsible for the rapid decrease of the latter term as the 
particle with increasing A moves in orbits with larger and 
larger radius. It should be emphasized that the energies 
(10) will be somewhat modified in a many mode excited 
state as a result of mode-mode interactions, but these 
interactions contribute terms that are at most of order 
v and therefore small compared with the leading order 
term (10). Thus, at least for L <C N, we can consider 
the Q\ excitations as a gas of non-interacting particles. 

Attractive interactions (v < 0) 

On the basis of the above, we are now in a position to 
construct the yrast states and the low energy excitations 
as a function of the angular momentum L. The phase dia- 
gram is somewhat simpler for attractive interaction than 
for repulsive and so we begin with attraction. The mode 
energies (10) are seen to be positive for all modes with 
the exception of A = 1 which has ei = 0; the vanishing 
of the mode energy for A = 1 can be simply understood 
from the fact that the operator Qi is proportional to the 
center of mass coordinate and thus does not affect the in- 
teractions of the particles which only depend on relative 
coordinates. Thus, in agreement with [3], we find that 
the yrast state is described by a wave function 

1-l~(Qi) L $o, (11) 

in which the excitation energy and angular momentum 
are entirely associated with the center of mass of the atom 
cloud moving around in the confining oscillator potential 
while the relative coordinates and interactions remain ex- 
actly as in the L = state. 

The state (11) can be seen to involve two distinct con- 
densates: partly the usual Bose- Einstein condensate that 
puts all the particles into the lowest state of the oscillator. 
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This condensate is, in the presence of a finite tempera- 
ture, melted by excitations that move particles to excited 
states in the oscillator and this results in a critical tran- 
sition at the Bose-Einstein condensation temperature [1] 



(12) 



In addition, the state (11) involves a condensate into a 
single mode (A = 1) out of the many degrees of freedom 
in partition space. The low energy degrees of freedom of 
the partition space are described by excitations involving 
the removal of particles from the A = 1 mode and corre- 
sponding excitation of the modes with A = 2, 3 4, ... L. 
Such excitations imply a melting of the partition space 
condensate and because of the (near) degeneracy of all 
the A ^ 1 normal modes, the thermal average of their 
contributions to the partition function can be be easily 
evaluated and yield a critical temperature 



Nv 

pc ~Jn~L 



(13) 



for the partition space phase transition. At this temper- 
ature the specific heat has a singularity (decreasing by 
a factor of 2) and the occupation of the (A = 1) mode 
decreases from the value 2L(1 + V5) 1 to a value that 
is of order yL, This transition is, of course, in a finite 
system and therefore the critical region is spread over a 
finite interval in temperature. However this interval is 



ST, 



Nv 



(14) 



(InL) 2 

and thus small compared to the critical temperature. 
Repulsive interactions (v > 0) 

The spectra resulting from repulsive interactions are ob- 
tained by inverting those discussed above and the struc- 
ture of the yrast region involves an analysis of the config- 
urations that would have been those of highest energy in 
the earlier discussion. It turns out that there are again 
two distinct condensates at each value of L and two phase 
transitions as a function of temperature, and in addition 
there is a third type of phase structure giving rise to a 
quantal phase transition (i.e. a transition at T = 0) in 
which the nature of the condensate in partition space 
changes abruptly as a function of L. These different 
T = phases are distinguished by the different integer 
number of vortex lines that exist in the ground state. 

We can again, for repulsive interaction (and L -C N), 
consider the spectra constructed from the collective ex- 
citations produced by the operators Q\, but now the ex- 
pression (10) yields negative values. The fact that the ex- 
citations have a negative energy reflects the fact that for 
repulsive interactions the interaction energy in the L = 
state is a maximum and all the Q\{\ ^ 1) excitations re- 
duce the energy by putting nodes in the wave function for 



relative motion and letting the particles get a little fur- 
ther away from each other. The yrast states will therefore 
be obtained by a condensation into that mode that has 
the greatest energy gain per unit of angular momentum. 
As can be seen from Table I, the quadrupole (A = 2) 
and octupole (A = 3) modes are optimal in this respect; 
thus, for L <C N, the yrast state will be a condensate 
involving only these two modes [8]. Excitations of the 
condensate are obtained by removing quanta from the 
condensate mode and placing them in other Q\ modes, 
always subject to the constraint (4). In the presence of 
a finite temperature such excitations deplete the conden- 
sate and lead to a phase transition at a temperature of 
order (13). 

As we go to higher angular momentum quite general 
arguments Q suggest that for repulsive interactions and 
L ~ N the near yrast states will be dominated by a 
structure involving a unit vortex line 



n 



N • 
p=l^P. 



$0- 



(15) 



In the state (15) all N atoms have been collectively 
shifted outward in the radial direction of the (x, y) plane 
effectively expanding the area, reducing the density, and 
thus reducing the energy associated with the short range 
repulsive interactions. 

The state (15) represents a realization of the partition 
(1^) in the mode corresponding to (5) and as discussed in 
connection with equation (7) represents a terribly com- 
plicated state if expressed in terms of Q\ excitations. 
Since there is no efficient way of combining a condensate 
in the Q2 mode with a condensation in the to = 1 single 
particle state (unit vortex) , the yrast state must undergo 
a violent rearrangement (many narrowly avoided cross- 
ings) at the critical angular momentum that separates 
these two regimes, ie a quantum phase transition. 

In the region around L = N the low energy degrees of 
freedom can be described by the Q\ operators acting on 
the vortex state, but now the excitation energies e\ are 
different from those give by (10) since the, A — excitation 
now represents a mode in which a particle is removed 
from the condensate in (to = 1) and excited to a state 
with to = A + l. There is also an additional mode in these 
spectra that results from the possibility of removing a 
particle from the to = 1 condensate and exciting it to 
the to = state. This A = — 1 mode can be excited by 
the operator 



1 N 1 

in at ^ — ' 7. 



(16) 



However the multiple excitations of this mode are not de- 
scribed by powers of the operator (16) since these pow- 
ers contain terms that go outside the partition space. 
Rather, the nth excitation of the A = — 1 mode is cre- 
ated by the operator 



E 

Pl<P2-<Pn 



(z pi Zp 



z p n y 



(17) 
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The degree of freedom associated with the A = — 1 mode 
can be seen as a displacement of the vortex line perpen- 
dicular to the axis along which the angular momentum is 
aligned as described by the broken symmetry condensate 
wave function 

* ~ (lTf =1 (z p - a)) *„• (18) 

The many interesting features of rotating Bose conden- 
sates that are briefly sketched in the present note, are the 
subject of a more systematic report that is under prepa- 
ration in collaboration with Christopher Pethick. In this 
connection it should be mentioned that the rather special 
features of the confining potential assumed in the present 
note, do not appear to be crucially important; rather the 
qualitative features of the phase diagram with its multi- 
plicity of phase transitions exhibit a considerable range 
of robustness under changes involving anharmonicity and 
anisotropy of the external potential. 

I am indebted to my colleague Christopher Pethick in 
more ways than one: partly for his inspiring introduction 
to the Wonderful World of Bose Condensates in external 
traps, partly for informing me of his results, obtained 
in collaboration with L.P. Pitaevskii, on the structure of 
the ground state for L ^ and attractive interaction, but 
most of all for his posing the seminal question, "what is 
the structure of the yrast line for these systems?" . I have 
also received valuable help from A.D. Jackson in address- 
ing the combinatorial problems encountered in construct- 
ing quantum states in partition space. I wish to thank 
Georgios Kavoulakis for pointing out an error in the orig- 
inal manuscript for this article and for his collaboration 
in the further development. 



[8] A detailed study, including the effect of mode-mode in- 
teractions in the low angular momentum region has been 
carried out in collaboration with G. Kavoulakis and C. 
Pethick. This work establishes the fact that in this region 
the condensate involves only the quadrupole mode until 
a critical angular momentum at which both octupole and 
quadrupole modes become involved. 
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£a/A 


1 








2 


-1/2 


-1/4 


3 


- 3/4 


-1/4 


4 


- 7/8 


- 7/32 


5 


- 15/16 


- 3/16 



Table 1: Caption: 

Energies of phonon modes for repulsive interactions. 
The energies (in units of -/Vjt)|) are obtained from (10), 
and apply to Q\ excitations in the yrast region for 
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